Dynamics of quantum phase transition: exact solution in quantum Ising model 
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Quantum Ising model is an exactly solvable model of quan- 
tum phase transition. This paper gives an exact solution when 
the system is driven through the critical point at finite rate. 
The evolution goes through a series of Landau-Zener level 
anticrossings when pairs of quasiparticles with opposite pseu- 
domomenta get excited with probability depending on the 
transition rate. Average density of defects excited in this way 
scales like a square root of the transition rate. This scaling is 
the same as the scaling obtained when the standard Kibble- 
Zurek mechanism of thermodynamic second order phase tran- 
sitions is applied to the quantum phase transition in the Ising 
model. 

PACS numbers: 03.65.-w, 73.43.Nq, 03.75.Lm, 
32.80.Bx, 05.70.Fh 

Introduction. — Phase transition is a fundamental 
change in a state of a system when one of parameters of 
the system passes through its critical point. The states 
on the opposite sides of the critical point are character- 
ized by different types of ordering. In a second order 
phase transition the fundamental change is continuous 
and the critical point is characterized by diverging cor- 
relation length and relaxation time. This critical slow- 
ing down implies that no matter how slowly a system 
is driven through the transition its evolution cannot be 
adiabatic close to the critical point. If it were adiabatic, 
then the system would continuously evolve between the 
two types of ordering. However, in the necessarily non- 
adiabatic evolution ordering of the state after the tran- 
sition is not perfect, the state is a mosaic of ordered do- 
mains whose finite size depends on the rate of the tran- 
sition. This scenario was first described by Kibble [1] 
and then its underlying dynamical mechanism was pro- 
posed by Zurek [1] who predicted that the size of the 
ordered domains scales with the transition time tq as 
Tq, where w is a combination of critical exponents. The 
Kibble-Zurek mechanism (KZM) of second order ther- 
modynamic phase transitions was confirmed by numeri- 
cal simulations of the time-dependent Ginzburg-Landau 
model [2] and tested by experimens in liquid crystals [3] , 
superfluid helium 3 [4] , superfluid helium 4 [5] , and both 
high-T c [6] and low-T c [7] superconductors. KZM is a 
universal theory of dynamics of second order phase tran- 
sition whose applications range from the low tempera- 
ture Bose-Einstein condensation (BEC) [8] to the ultra 
high temperature transitions in the grand unified theo- 
ries of high energy physics. However, the zero tempera- 
ture quantum limit remains largely unexplored (but see 
Refs. [9-11]) and quantum phase transitions are in many 
respects qualitatively different from transitions at finite 



temperature. Most importantly time evolution is unitary, 
there is no damping that seems to play essential role in 
KZM. 

According to Sachdev [12] understanding of quantum 
phase transitions is based on two prototypical models. 
One is Bose- Hubbard model with its transition from Mott 
insulator to superfluid and the other is the 1-dimensional 
quantum Ising model. Of the two only the Ising model 
is exactly solvable. It is defined by the Hamiltonian 



H 



N 



n u n+l) 



with periodic boundary conditions 



(1) 



(2) 



The 1-dimensional quantum Ising model has the same 
critical properties as the 2-dimensional classical Ising 
model. Quantum phase transition takes place at the crit- 
ical value g = 1 of external magnetic field. When g ^> 1, 
then the ground state is a paramagnet | —>■—>— ► . . . — ►) 
with all spins polarized up along the x-axis. On the other 
hand when g <C 1, then there are two degenerate ferro- 
magnetic ground states with all spins pointing either up 
or down along the z-axis: | f"TT • • • T) or I III •••!)• I n 
adiabatic transition from paramagnet to ferromagnet the 
system would choose one of the two ferromagnetic states 
(with possible help of infinitezimal symmetry breaking 
field along the z-axis, as usual in symmetry breaking 
phase transitions). However, when N — > oo, then en- 
ergy gap at g = 1 tends to zero (quantum version of the 
critical slowing down) and it is impossible to pass the 
critical point without exciting the system. As a result 
the system ends in a quantum superposition of states 
like 



Tiiiiimmniiimmi. 



(3) 



with finite domains of spins pointing up or down and sep- 
arated by kinks where the polarization of spins changes 
its orientation. Average size of the domains or, equiv- 
alently, average density of kinks depends on transition 
rate. When transition is slow, then the size is large, 
but when it is very fast, then orientation of individual 
spins can become random, uncorrelated with their near- 
est neighbors. Transition time tq can be unambiguously 
defined when we assume that close to the critical point 
at g = 1 time-dependent field g(t) driving the transition 
can be approximated by a linear quench 



1 



g(t<0) = - 



(4) 



with variable quench time tq . Density of kinks after the 
linear quench was estimated in Rcf. [11] as 



2Jt q 



1/2 



(5) 



in an attempt to generalize KZM to quantum phase tran- 
sitions. This is an order of magnitude estimate with an 
unknown 0(1) prefactor. Results of numerical simula- 
tions in Ref. [11] for a chain of N spins were fitted with 
n <~ Tq W . The fitted exponent w decreased with increas- 
ing N reaching w = 0.58 for N = 100. This is roughly 
consistent with the predicted w = \. 

These encouriging results strongly motivate further 
study of the dynamics of quantum phase transition in the 
quantum Ising model. As this model is exactly solvable 
dynamics of its quantum phase transition also deserves 
to be solved exactly. In fact this is the only model, with- 
out any counterpart in the thermodynamic transitions, 
where KZM can be compared with an exact analytic so- 
lution. I begin by reminding relevant facts about energy 
spectrum of the model. 

Energy spectrum. — Here I assume N is even for 
convenience. After Jordan- Wigner transformation [13], 



(6) 
(7) 



introducing fermionic operators c„ which satisfy anti- 
commutation relations |c m ,c^} = 5 mn and {c m ,c„} = 
{cj^c^} = the Hamiltonian (1) becomes [14] 



H 



p + H + p+ + p- H - p- 



where 



N 



i ± n < 



i ± n (i - 24 c) 



(9) 



even parity. H + is diagonalized by Fourier transform fol- 
lowed by Bogoliubov transformation [14]. Fourier trans- 
form consistent with the antiperiodic boundary condition 
CN+l = -ci 



3 — iir/4 

~7W 



E 



c k e 



ik(na) 



(ii) 



where the pseudomomenta k take "half integer" values 

(12) 



k ~ ± 2Na' 



,±1^-1 
2 2 



2tt 



transforms the Hamiltonian into 

H+ = 

^2[g — cos(ka)]c\ck + sin(fca) 



C k C -k + C -k c k 



-g} ■ (13) 



Here a is lattice spacing. Diagonalization of H + is com- 
pleted by the Bogoliubov transformation 



Cfc 



(14) 



provided that Bogoliubov modes (uk,Vk) are eigenstates 
of the stationary Bogoliubov-de Gennes equations 



e Uk = +2J[.g — cos(ka)]uk + 2Jsin(fca)ffe 
e Vk — —2J[g — cos(ka)]vk + 2 J sin(ka)uk 



(15) 



There are two eigestates for each k with cigenencrgies 
e = ±efe, where 



e k = 2J^[g - cos(ka)} 2 + sin 2 (ka) . (16) 
The positive energy eigenstate 



(8) (uk, Vk) ~ (g — cos ka) + \/g 2 — 2g cos ka+1, sin(fca) 



normalized so that |wfc| 2 + \vk\ 2 = 
quasiparticle operator -f k = u k c k 



(17) 

1 , defines the 



V-kC 



-k ' 



and 



the negative energy eigenstate (u k ,v k ) = (—Vk,Uk) 



are projectors on the subspaces with even (+) and odd 
(— ) numbers of c-quasiparticles and 

N 

ff ± = J ^2 (gcicn - 4c„+i - c„+ic„ - | + h-c) ■ (10) 

n=l 

are corresponding reduced Hamiltonians. The c„'s in H~ 
satisfy periodic boundary conditions c^+i = c\, but the 
c„'s in H + must obey cjy+i = — C\ what I call "antiperi- 
odic" boudary conditions. 

The parity of the number of c-quasiparticles is a good 
quantum number and the ground state has even par- 
ity for any value of g. Assuming that a quench begins 
in the ground state we can confine to the subspace of 



defines lk = {u k )*c k + v_ k d k = - 7 I fc . After 
the Bogoliubov transformation the Hamiltonian H + = 

\ Efc £ fe (ihk - Ikhk) equivalent to 



H+ 



E £fc 

k 



7fe7fc 



(18) 



which is a simple-looking sum of quasiparticles with half- 
integer pseudomomenta. However, thanks to the projec- 
tion P + H + P + in Eq.(8) only states with even numbers 
of quasiparticles belong to the spectrum of H. 

Transition from paramagnet to ferromagnet. — 
In the linear quench Eq.(4) the system is initially pre- 
pared in its ground state at large initial value of g >• 1 



2 



but when g is ramped down to zero the state of the sys- 
tem gets excited from its instantaneous ground state and 
in general its final state at t = has finite number of 
kinks. Comparing the Ising Hamiltonian Eq.(l) at g = 
with the Bogoliubov Hamiltonian (18) at g = we ob- 
tain a simple expression for the operator of the number 
of kinks 



TV 



1 N 



k 



(19) 



The number of kinks is equal to the number of quasipar- 
ticles excited at g — 0. The excitation probability 



(20) 



in the final state can be found with the time-dependent 
Bogoliubov method. 

The initial ground state is Bogoliubov vacuum |0) an- 
nihilated by all quasiparticle operators j k which are de- 
termined by the asymptotic form of the (positive en- 
ergy) Bogoliubov modes (u kl v k ) ~ (1,0) in the regime 
of g 3> 1. When g{t) is ramped down, then the quantum 
state | ip(t)} in general gets excited from the instantaneous 
ground state. The time-dependent Bogoliubov method 
makes an Ansatz that \ip{t)} is Bogoliubov vacuum anni- 
hilated by a set of quasiparticle annihilation operators ^ k 
defined by a time-dependent Bogoliubov transformation 



Cfe 



u k (t)% +v*_ k (t)^_ k 



(21) 



with the initial condition [uk{— 00), v k (— 00)] = (1,0). 
The Bogoliubov modes [u k (t), Vk(t)] must satisfy Heisen- 
berg equation ih-j^c k — [c k ,H + ] with the constraint 

3j7fc = [equivalent to 7fc|"0W) = 0]- The Heisenberg 
equation is equivalent to the dynamical version of the 
Bogoliubov-de Gennes equations (15): 



ifr—Uk = +2 J [g(t) — cos(fca)] u k + 2Jsin(fca) v k , 

zft^Vfc = — 2J [g(t) — cos(fca)] v k + 2Jsin(/ca) u k ■ (22) 

At any value of g Eqs.(22) have two instantaneous eigen- 
states. Initially the mode [u k (t), v k {t)\ is the positive en- 
ergy eigenstate, but during the quench it gets "excited" 
to a combination of the positive and negative mode. At 
the end of the quench at t = when g = we have 



Pk = \Pk\ 



(25) 



which can be found by mapping Eqs.(22) to the Landau- 
Zener (LZ) problem (similarity between KZM and LZ 
problem was first pointed out by Damski in Ref. [16]). 
The transformation 

t = 4Jtq sin(fca) ( — + cos(fca) ) (26) 
\ T Q ) 

brings Eqs.(22) to the standard LZ form [15] 



■ft d 
in—Uk 

dr 



1 , . » 1 
■-(rAfe)Mfe + -Vk 



d 1 - . . 1 



(27) 



with A^ 1 = 4 Jtq sin 2 (ka). Here the time r runs from 
—00 to Tfinai = 2 Jtq sin(2fca) corresponding to t = 0. 
Tunnelling between the positive and negative energy 
eigenstates happens when r e (— A^ 1 , A^ 1 ). Tfi na i is 
well outside this interval, Tfinai 3> A^T , for long wave- 
length modes with \ka\ <C f . For these modes time r 
in Eqs.(27) can be extended to +00 making them fully 
equivalent to LZ equations [15]. 

In slow transitions only long wavelength modes can get 
excited. For these modes we can use the LZ formula [15] 
for excitation probability: 



Pk 



-2tt( Jtq /h)(ka) 2 



(28) 



This calculation is self-consistent when the width of the 
obtained gaussian ka = (^ttJtq/H)^ 1 ^ 2 is much less than 
j or, equivalently, for slow enough quenches with tq ^> 
■^rj. With the LZ formula (28) we can calculate the 
number of kinks in Eq.(19) as 



Pk 



(29) 



There are at least two interesting special cases. 

1) When N — > 00 i.e. in the limit of proper phase tran- 
sition the sum in Eq.(29) can be replaced by an integral. 
Expectation value of density of kinks becomes 



n = hm — 



1 

27 



d(ka) p k 



— 1/2 

The density scales like Tq in agreement with KZM. It 



1 



1 



2tt x /2Jt q /H' 



(30) 



K(0), v k (0)} - a k (u k ,v k ) + [3 k [u k , v k ) (23) ig & factor f 1 ~ .159 less than the simple estimate 



and consequently 7^ = a k ^ k — Pkl k - The final state which 
is, by construction, annihilated by both j k and j- k is 

W0)) = n(«*+^7h-fc)|0>- (24) 

fc>0 

Pairs of quasiparticles with pseudomomenta (k,—k) are 
excited with probability 



in Eq.(5) based on KZM [11]. The results of numerical 
simulations in Ref. [11] are best fitted with a prefactor of 
0.16 which is consistent with the exact number. 

2) For a finite chain we can ask what is the fastest tq 
when still no kinks get excited. This critical Tg d marks a 
crossover between adiabatic and non-adiabatic regimes. 
In other words we can ask what is the probability for 
a finite chain to stay in the ground state. As different 



3 



pairs of quasiparticlcs (k, —k) evolve independently, the 
probability to stay in the ground state is the product 

fc>0 

Well on the adiabatic side only the pair (77,-77) is likely 
to get excited and we can approximate 

V GS « 1- Pj r « l-exp^Tr 3 ^ . (32) 

Again, numerical results in Ref. [11] are consistent with 
this formula: the best fit gives 59 in place of the 2tt 3 62. 
A quench in a finite chain becomes non-adiabatic when 
tq is less than Tg d = t^wj which grows with the system 
size like N 2 . In other words, the size N of a defect-free 

1 /2 

chain grows like Tq in consistency with Eq.(30). 

Transition from ferromagnet to paramagnet. — 

Again, close to the critical point at g = 1 the ex- 
ternal field can be approximated by a linear quench 
git > 0) = The initial state is the even parity ground 
state at g = and the final state \i>(i)} is in general 
excited with respect to the polarized ground state at a 
large final j> 1. We want to know how many spins in 
the final state are flipped with respect to the polariza- 
tion in the ground state or, more precisely, what is the 
expectation value of the operator of the number of flips 

T = \ (N - J2n=i a n) ■ Wc notc that whcn 9 > !> thcn 

H « ~JgT,n=i< = ~ J 9( N - 2-T 7 ), but on thc othcr 
hand H + w 2JgJ2k7k7 k ~~ ^9^, and consequently 

T = £^fe ■ (33) 
k 

At large j> 1 the number of flips is the number of ex- 
cited quasiparticles. Formal similarity of Eqs.(33) and 
(19) allows us to proceed in the same way as in thc para- 
magnet to ferromagnet transition. We soon arrive at 
Eqs.(27) but with r replaced by — t on thc right hand 
side. In analogy to Eq.(30) we find density of flipped 
spins 

1 = 2^ ^2Jr Q /h (M) 

and the probability for a finite chain to stay in the ground 
state is again given by Eq.(31). 

Conclusion. — This paper gives an exact analytic so- 
lution when the quantum Ising model is driven through 
its quantum critical point at finite transition rate. The 
evolution goes through a series of Landau-Zener level 
anticrossings when pairs of quasiparticles with opposite 
pseudomomenta get excited with probability depending 
on the transition rate. Average density of defects excited 
in this way scales like a square root of the transition 



rate. This scaling is the same as the scaling obtained 
when the standard Kibble-Zurek mechanism of thermo- 
dynamic second order phase transitions is extended to 
the quantum phase transition in the Ising model. 
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